Thermal conductance of graphene and dimerite 
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We investigate the phonon thermal conductance of graphene regarding the graphene sheet as the 
large- width limit of graphene strips in the ballistic limit. We find that the thermal conductance 
depends weakly on the direction angle 9 of the thermal flux periodically with period 7r/3. It is 
further shown that the nature of this directional dependence is the directional dependence of group 
velocities of the phonon modes in the graphene, originating from the D^h symmetry in the honey- 
comb structure. By breaking the D^h symmetry in graphene, we see more obvious anisotropic effect 
in the thermal conductance as demonstrated by dimerite. 

PACS numbers: 81.05.Uw, 65.80.-|-n 



I. INTRODUCTION 

As a promising candidate material for nanoelectronic 
device, graphene has been attracted intensive atten- 
tion in research in past years (for review, see e.g. 
Ref. [H). It demonstrates not only peculiar electronic 
properties^i^, but also very high (as high as 5000 
Wm~-'^K~^) thermal conductivity,^'^ which is beneficial 
for the possible electronic and thermal device applica- 
tions of graphene i^i^i^ii^ A recent work has studied the 
structure of an interesting new allotrope of graphene by 
a first principle calculation.^^ The ground state energy 
in this allotrope is about 0.28 eV/atom above graphene, 
which is 0.11 eV/atom lower than Ceo- 

In this paper, we calculate the ballistic phonon ther- 
mal conductance for the graphene sheet by treating the 
graphene as the large width limit of graphene strips, 
which can be described by a lattice vector R — nidi + 
712(22 »^ The phonon dispersion of the graphene is ob- 
tained in the valence force field model (VFFM), where 
the out-of-plane acoustic phonon mode is a flexure mode, 
i.e., it has the quadratic dispersion around F point in 
the Brillouin zoneJ^ Our result shows that the thermal 
conductance has a T^-^ dependence at low temperature, 
which is due to the contribution of the flexure modeJ^ 
At room temperature, our result are comparable with the 
recent experimental measured thermal conductivityi^ 

We find that the thermal conductance in graphene de- 
pends on the direction angle 9 of the thermal flux pe- 
riodically with 7r/3 as the period. The difference be- 
tween maximum and minimum thermal conductance at 
100 K is 1.24xl0^Wm-2K-i, which is about 1% varia- 
tion. Our study shows that this directional dependence 
for the graphene is attributed to the directional depen- 
dence of the velocities of the phonon modes, which origins 
from the Dgh symmetry of the honeycomb structure. 

For the dimerite, where the ZJgft symmetry is broken, 
the thermal conductance shows more obvious anisotropy 
of 10% and the value is about 40 % smaller than that of 
the graphene at room temperature. 



The present paper is organized as follows. In Sec. II, 
we describe the graphene strip by a lattice vector. The 
formulas we used in the calculation of the thermal con- 
ductance are derived in Sec. III. Calculation results for 
graphene and dimerite are discussed in Sec. IV A and 
Sec. IV B, respectively. Sec. V is the conclusion. 

II. CONFIGURATION 

In graphene, the primitive lattice vectors are ai and 
a2, with |ai| — |a2| — V3&o- — 1-42 is the C-C bond 
length in graphene. The corresponding reciprocal unit 
vectors are bi ^ (|l, ^ _^). 

As shown in Fig.[T](a), a strip in the graphene sheet can 
be described by a lattice vector R = nidi+n2d2. The real 
lattice vector H = pidi + P2d2 is introduced through^^: 
nip2 — n2Pi = N (N is the greatest common divisor of 
ni and 712)- The strip is denoted by NhH x NrR, where 
Nh and Nr are numbers of the periods in the directions 
along H and R, respectively. Instead of oi and a2, we use 
(i/, R/N) as the basic vectors in the following, and bn 
and bfi are their corresponding reciprocal unit vectors: 

bH = j^{-n2bi + riih), 

Any wave vector in the reciprocal space can be written 
as: 

k = knbH + kubji. (1) 

Using the periodic boundary conditions, this strip has 
Nh translational periods in the H direction and Nr x N 
translational periods in the R direction. As shown in 
Fig. [T] (b), the Brillouin zone for the graphene strip is 
Nh X N discrete segments, which are parallel or co- 
incide with bu- The coordinates for the wave vectors 
on these lines arei^ (fc/f, kii)~{i/NH, j/N^N), with 
i = 0, 1, 2, Nh - 1 and j = 0, 1, 2, ...,NbN - 1. 
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FIG. 1: Graphene strip is described by a lattice vector R = 
nioi + n2a2. (a). Strip with (ni,n2) = (4,2) and {Nh, 
Nr)={12, 1); (b). The Brillouin zone for this special strip is 
two discrete segments (solid) in the reciprocal space. 



(2) 



where 9 determines the direction of the thermal flux: 

-hR is the wave vector 



NNi 



eg = (cos 6*, sin 6'). k = kubH 
in the Brillouin zone of the strip, i.e., on the Nr. x N dis- 
crete lines. The transmission function for a phonon mode 

Tn{k) is assumed to be one."'^^ '^n(fc) = '^^dke'' ^"'^^ group 
velocity of mode (fc, n) in eg direction. The value of the 
group velocity can be accurately calculated through the 
frequency and the eigen vector of this phonon mode i^^i^° 
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where D is the dynamical matrix and Un{k) is the eigen 
vector. Only those phonon modes with v^> Q contribute 
to the thermal conductance in the eg direction. 

In the two-dimensional graphene strip system, it is 
convenient to use conductance reduced by cross section: 
a = a/ where s — W x h \s the cross section. The 
thickness of the strip, h — 3.35 A, is chosen arbitrarily 
to be the same as the space between two adjacent layers 
in the graphite. The width for the strip is W = Nii\R\, 
where, the thermal flux in the strip is set to be in the 
direction perpendicular to R, i.e., eg = bnlbn- We ad- 
dress a fact that the integral parameter (kn ) in Eq. (0) 
is the quantum number along the thermal flux direction. 

The thermal conductance in eg direction of the 
graphene can be obtained by: 



ct(T) = lim — cr(r). 



(4) 



The graphene sheet is actually a strip in the limit of 
Nh — > oo and — > oo. In this case, the Brillouin 
zone for the strip, i.e., x N discrete lines, turns to 
the two-dimensional Brillouin zone for the graphene. 

III. CONDUCTANCE FORMULAS 

The contribution of the phonon to the thermal con- 
ductance in the ballistic region is; i6,i7^i8 

1 /"^ df 
a{T) = ^ T{uj)huj—duj, 

where f{T,u}) is the Bose-Einstein distribution function. 
T{uj) is transmission function. In the ballistic region, 
T{u)) is simply the number of phonon branches at fre- 
quency LO. 

From the above expression, the thermal conductance 
in the graphene strip can be written as: 

NNn-l 6 .bH 
j=0 n=lvi)>0 " 



IV. CALCULATION RESULTS AND 
DISCUSSION 

A. graphene results 

The phonon spectrum of graphene is calculated in 
the VFFM, which has been successfully applied to 
study the phonon spectrum in the single-walled carbon 
nanotubesi^ and multi-layered graphene systemsi^ In 
present calculation, we utilize three vibrational poten- 
tial energy terms. They are the in-plane bond stretch- 
ing (V;) and bond bending (Vbb), and the out-of-plane 
bond bending (K-c) vibrational potential energy. The 
three force constants are taken from Ref. [2ll as: fc;=305.0 
Nm-\ /cbb=65.3 Nm"! and fcBB = 14.8 Nm"!. 

1. temperature dependence for thermal conductance 

In Fig. [21 the temperature is 100 K and the direction 
angle for the thermal flux is 6* = 7r/3. It is shown that the 
thermal conductance for a strip decreases with increasing 
width. At about VF=100 A, the thermal conductance 
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TABLE L The dependence of the thermal conductance on group velocities of phonon modes, 'b' in the 2nd line is a fitting 
parameter (see text). In the 4th line, the down (up) arrow indicates the decreasing (increasing) of Act when the corresponding 
phonon mode is excited more. 



velocity 


a 


V2 


V3 


V4 


V5 


^6 


sign(6) 






+ 




+ 




5" (X 


1 


1 


1 

"a 








AS 


I 


i 


T 


T 


I 


T 



CM 



03 



1.5 r 

1.4 - 
1.3 
1.2 
1.1 
1 



T=100 K 







100 200 
W (A) 



300 



400 



FIG. 2: Convergence for the thermal conductance of the 
graphene strip with increasing width at temperature 100 K. 
In the large width limit (VF>150A) the thermal conductance 
of the strip can be considered as the thermal conductance 
value of graphene. 
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FIG. 3: The thermal conductance of the graphene sheet v.s. 
temperature. Inset is logo- v.s. logT in extremely low tem- 
perature region. The calculated results (filled squares) can be 
fitted by function f{x) = 13.44 + 1.5a; (blue line). It indicates 
that the thermal conductance has a T^'^ dependence in this 
region. 



reaches a saturate value, which is actually the thermal 
conductance for the graphene. In the calculation, the 
width we used is about 300 A, which ensures that the 
strip is wide enough to be considered as a graphene sheet. 

In Fig. [31 the thermal conductance versus the temper- 
ature is displayed. In the low temperature region, the 
thermal conductance has a T^-^ dependence. This is the 
result of the flexure mode in the graphene sheet, which 
has the dispersion u = ak^. In the very low tempera- 
ture region, this mode makes the largest contribution to 
the thermal conductance. Its contribution to the thermal 
conductance is^^ a oc T^'^/^/a, which can be seen from 
the figure in the low temperature region. At room tem- 
perature T=300 K, the value for the thermal conductance 
is about 4.4 x 10^ Wm~^K~-'^. This result agrees with the 
recent experimental value for the thermal conductance in 
the graphene4i^ In the experiment, the thermal conduc- 
tivity is measured to be about 5.0 x 10^ Wm~^K~-'^ at 
room temperature. The distance for the thermal flux 
to transport in the experiment is L=11.5 ^m. So the 
reduced thermal conductance can be deduced from this 
experiment as tr = f = |; = 0.43 x 10^ Wm~'^K~^. Our 
theoretical result is much larger than this experimental 
value. Because our calculation is in the ballistic region, 
while in the experiment, there is scattering on defects, 
edges or impurities and thus the transport is partially 
diffusive. At the high temperature limit T = lOOOiir, our 
calculation gives the value 8.9 xlO^ Wm^^K^^, which is 
in consistency with the previous theoretical result 



2. directional dependence for thermal conductance 

As shown in Fig. [U at r=lGO K, the thermal conduc- 
tance varies periodically with the direction angle 9. The 
calculated results can be fitted very well by the function 
f{9) = a + bcos{66) -I- ccos(126'). The difference between 
the thermal conductance in the two directions with an- 
gle 6* = and n/2 is about 1.2x lO'^Wm-^K-i. This 
difference is very stable for graphene strips with differ- 
ent width (see Fig. E]). At T=100 K, the lattice thermal 
conductance is about two orders larger than the elec- 
tron thermal conductancei^ So the experimental mea- 
sured thermal conductance at T=100 K is mainly due to 
the contribution of the phonons. As a result, our calcu- 
lated directional dependence of the lattice thermal con- 
ductance in the graphene can be carefully investigated in 
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FIG. 4: The direction dependence of thermal conductance. 
6 is the direction angle for the thermal flux. The calculated 
results (filled squares) are fitted by the function f(x) = a + 
6cos(66i) + ccos(126') with a = 1.0456 x 10^ h = -6.237 x 10^ 
and c = -9.642 x 10^ 
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FIG. 5: Thermal conductance in 6 = Q and 7r/2 directions 
versus width of graphene strip. Lines are guide to the eye. 



the experiment. In the following, we say that two quan- 
tities Qi = Oi + hi cos Q9 and Q2 — a2 + ^2 cos &9 have 
the same (opposite) dependence on 6*, if the signs of bi 
and 62 are the same (opposite). 

To find the underlying mechanism for this directional 
dependence for the thermal conductance, firstly we show 
in Fig. [S] the coefficient a for the fiexure mode and the 
velocities for the other five phonon modes at the T point. 
Interestingly, this coefficient and velocities are also direc- 
tional dependent with the period 7r/3. Obviously, they 
can be fitted by function f{9) — a + hcos{(!>9). In Tabled 
the sign of the fitting parameter h for this coefficient and 
five velocities are listed, which can be read from Fig. [6] 
In the third line of Table U we list the contribution of 
the six phonon modes to the thermal conductance. If 
the three low frequency modes are excited, the thermal 



FIG. 6: The coefficient and velocity of the six phonon spec- 
trum around T point in the Brillouin zone: (a) . the coefficient 
of the out-of-plane acoustic mode with uj — ak^ in the unit 
of 10~® m^s~^; (b)-(f). velocities of the other five phonon 
spectrum (from low frequency to high frequency), in the unit 
of ms~^ . The horizontal axises in all figures are the direction 
angle 9 in the unit of n. 
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FIG. 7: The difference of the thermal conductance between 
directions of 6 — and 7r/2 versus temperature. This quan- 
tity has an abundant temperature dependence. Inset is the 
enlarged figure for the low temperature region. 



conductance is in inverse proportion to their velocities 
While as can be seen from Eq. ([4]), the thermal conduc- 
tance is proportional to the velocities for the three high 
frequency optically modes when they are excited. In each 
temperature region, there will be a key mode which is the 
most important contributor to the thermal conductance. 
The direction dependence of the velocity of this key mode 
determines the direction dependence of the thermal con- 
ductance. 

We then further study the difference between the ther- 
mal conductance in two directions with 9 — and ^: 
Act = (T(f ) - ct(0). In the fourth fine of TablelH we dis- 
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play the effect of different modes on Act. It sliows tliat 
Act will decrease, if the first, second and fifth phonon 
modes are excited sufficiently with increasing tempera- 
ture. The other three phonon modes have the opposite 
effect on the thermal conductance. The dependence of 
Act on the temperature is shown in Fig. [71 where five 
different temperature regions are exhibited. 

(1) [0, 4]K: In this extremely low temperature region, 
only the flexure mode is excited. This mode results in 
Act < 0. Because the coefficient a depends on the di- 
rection angle 9 very slightly, the absolute value of Act is 
pretty small (see inset of Fig. [7]). 

(2) [4, 10]K: The second acoustic mode is excited in 
this temperature region. In respect that this mode has 
more sensitive direction dependence and favors to de- 
crease Act, Act decreases much faster than region (1). 

(3) [10, 70]K: In this temperature region, the third 
acoustic mode begins to have an effect on the thermal 
conductance. This mode's directional dependence is op- 
posite of the previous two acoustic modes and it will 
increase Act. The competition between this mode and 
the other two acoustic modes slow down the decrease of 
the value Act at temperature below T—4QK. The third 
acoustic mode becomes more and more important with 
temperature increasing, and Act begins to increase after 
T=40K as can be seen from the inset of Fig. [T] 

(4) [70, 500]K: The third acoustic mode becomes the 
key mode in this temperature region. As a result, Act 
changes into a positive value and keeps increasing. 

(5) [500, lOOOjK: In this high temperature region, the 
optical mode will also be excited one by one in the fre- 
quency order with increasing temperature. Since there 
are two optical modes (1st and 3rd optical modes) favors 
to increase Act, while only one optical mode (2nd opti- 
cal mode) try to decrease Act, the competition result is 
increasing of Act in the high temperature region. 

r=100K is in region (4), where the direction depen- 
dence of the thermal conductance is controlled by the 
velocity of the third mode, so the dependence of ct on 
in Fig. |4]is opposite to the dependence of velocity ^3(6') 
in Fig.iKc). 



B. dimerite results 

The adatom defect is used as a basic block to man- 
ufacture a new carbon allotrope of graphene, named 
dimeriter^ and the relaxed configuration of this new ma- 
terial is investigated by a first principle calculation. The 
unit cell for the dimerite is shown in Fig.[S](a), where the 
two basic unit vectors are ai and 02, with |ai| = 7.60 A, 
and |a2| = 7.05 A. The angle between these two vectors is 
0.77r. In each unit cell, there is a (7-5-5-7) defect, which 
leads to two anisotropic directions: the (7-7) direction 
and the (5-5) direction. (7-7) direction is from the center 
of a heptagon to the opposite heptagon, and the (5-5) 
direction is from the center of a pentagon to the neigh- 
boring pentagon. These two directions are perpendicular 
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FIG. 8: (Color online) Configuration for the dimerite. (a) is 
one unit cell for dimerite. There are 16 atoms in each unit cell. 
The length for the two unit vectors are 7.60 Aand 7.05 A. The 
angle between them is 0.7 vr. (b) shows 49 unit cells together. 
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FIG. 9: (Color online) Phonon spectrum for the dimerite 
along high symmetry lines in the Brillouin zone. Around F 
point, there are low frequency optical modes with frequency 
about 200 cm~^, which do not exist in pure graphene. 

to each other, and corresponding to 9 — QAtt, O.Qtt re- 
spectively in Fig. [5] (a), where is the direction angle 
with respect to x axis. From graphene to dimerite, the 
symmetry is reduced from Dgh to D2h- 

1. phonon dispersions in dimerite 

The phonon dispersion for the dimerite from the above 
VFFM is shown in Fig. [5) Similar with pure graphene, 
there are three modes with zero frequency. Two of them 
are the acoustic modes in the xy plane. The other zero- 
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FIG. 10: (Color online) The thermal conductance for the 
dimerite (blue solid line) and graphene (red dotted line) are 
compared in a large temperature region. The value of the 
dimerite is about 40% smaller than the graphene at room 
temperature. 
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FIG. 11: (Color online) The directional dependence for the 
thermal conductance of dimerite. The calculated results 
(filled squares) are fitted by function f{x) = a + bcos{2{9 + c)) 
with a=1.9xl0^, 6=1.1x10®, and c= O.Itt (blue line). 

frequency mode is a flexure mode with parabolic dis- 
persion of u! = l3k'^. This flexure mode corresponds to 
the vibration in the z direction. Different from the pure 
graphene, there are a lot of optical phonon modes with 
frequency around 150 cm^^ in the dimerite. 

2. thermal conductance in dimerite 

In Fig. [ini we compare the thermal conductance of the 
dimerite with the graphene. At room temperature, the 
thermal conductance of the dimerite is about 40% smaller 
than the graphene. As can be seen from Fig. [TTl the 
thermal conductance in the dimerite is more anisotropic 



than graphene. The calculated results can be best fitted 
with function f{x) = a + b cos(2(6' + c)) with a=1.9x 10^, 
b=l.lxlO^, and c= O.Itt. The thermal conductance has 
a minimum value at 9 — OAtt, and maximum value at 
9 = O.Ovr. As mentioned previously, the directions with 
9 = OAtt, O.Ovr are the (7-7) and the (5-5) directions, 
respectively. So, the thermal conductance in (5-5) direc- 
tion is 12% larger than (7-7) direction, which is about 
one order larger than that of the pristine graphene. We 
expect this anisotropic effect detectable experimentally if 
the precision in the current experiment can be improved. 



V. CONCLUSION 

In conclusion, we have calculated the phonon thermal 
conductance for graphene in the ballistic region, by con- 
sidering the graphene as the large width limit of graphene 
strips. The calculated value for the thermal conductance 
at room temperature is comparable with the recent ex- 
perimental results, while at high temperature region our 
results are consistent with the previous theoretical cal- 
culations. We have found that the thermal conductance 
is directionally dependent and the reason is the direc- 
tional dependence of the velocities of different phonon 
modes, which can be excited in the frequency order with 
increasing temperature. By breaking the Deh symmetry 
in graphene, we can see more obvious anisotropic effect 
of the thermal conductance as demonstrated by dimerite. 

We have following two further remarks: 

(1) . Since the anisotropic effect of thermal conduc- 
tance in graphene is small (1%), it requires high accu- 
racy in the calculation of phonon mode's group velocity 
to see this anisotropic effect. Thanks to the superiority 
of the VFFM, we can derive an analytic expression for 
the dynamical matrix and calculate accurately the value 
of the group velocity following Eq. ([3]). Thus we can ob- 
tain the 1% anisotropic effect in graphene as discussed in 
this manuscript. We also used the Brenner empirical po- 
tential implemented in the "General Utility Lattice Pro- 
gram" (GULP)22 to calculate the phonon dispersion and 
group velocity in graphene. For lack of analytic expres- 
sion for the dynamical matrix, we find that the accuracy 
is not high enough for the group velocity and it is diffi- 
culty to see this anisotropic effect. 

(2) . From symmetry analysis^'*, when thermal trans- 
port is in the diffusive region where the Fourier's law 
exists, the Dgh symmetry of graphene constrains the 
thermal conductivity to be a constant value. So the 
anisotropic effect in graphene can not be seen if the 
thermal transport is in the diffusive region, yet it can 
only be seen in the ballistic region as discussed in this 
manuscript. But the situation changes in dimerite, where 
the Deh symmetry is broken into D2h- The D2/1 symme- 
try does not constrain the thermal conductivity to be a 
constant value even the Fourier's law is valid. So we can 
expect to see the isotropic effect in the dimerite both in 
the ballistic and diffusive region. 
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Note added in proof. We have learned about a re- 



cent theoretical study on the thermal conduction in sin- 
gle layer graphene. Nika et al^ investigated the effects 
of Umklapp, defects and edges scattering on the phonon 
thermal conduction in graphene. Our result for the ther- 
mal conductivity is larger than their value. This is be- 
cause our study is in the ballistic region, which gives an 
upper limit for their calculation results. 
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